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ABSTRACT 

The rapid circularization and synchronization of the stellar components in an eccentric binary 
system at the onset of Roche lobe overflow is a fundamental assumption common to all binary stellar 
evolution and population synthesis codes, even though the validity of this assumption is questionable 
both theoretically and observationally. Here we calculate the evolution of the orbital elements of 
an eccentric binary through the direct three-body integration of a massive particle ejected through 
the inner Lagrangian point of the donor star at pcriastron. The trajectory of this particle leads to 
three possible outcomes: direct accretion onto the companion star within a single orbit, self-accretion 
back onto the donor star within a single orbit, or a quasi-periodic orbit around the companion star, 
possibly leading to the formation of a disk. We calculate the secular evolution of the binary orbit in 
the first two cases and conclude that direct impact accretion can increase as well as decrease the orbital 
semi-major axis and eccentricity, while self-accretion always decreases the orbital semi-major axis and 
eccentricity. In cases where mass overflow contributes to circularizing the orbit, circularization can 
set in on timescales as short as a few per cent of the mass transfer timescale. In cases where mass 
overflow increases the eccentricity, the orbital evolution is governed by competition between mass 
overflow and tidal torques. In the absence of tidal torques, mass overflow results in direct impact can 
lead to substantially subsynchronously rotating donor stars. Contrary to assumptions common in the 
literature, direct impact accretion furthermore docs not always provide a strong sink of orbital angular 
momentum in close mass-transferring binaries; in fact we instead flnd that a significant part can be 
returned to the orbit during the particle orbit. The formulation presented in this paper together 
with our previous work can be combined with stellar and binary evolution codes to generate a better 
picture of the evolution of eccentric, Roche lobe overflowing binary star systems. 
Subject headings: Celestial mechanics. Stars: Binaries: Close, Stars: Mass Loss, Accretion, Methods: 
Numerical 



1. INTRODUCTION 

In close binaries, the combined effects of stellar and 
orbital evolution can cause a star to fill its Roche lobe 
and transfer mass to its companion. If the orbit is ec- 
centric this mass transfer is expected to commence at 
or near the periastron of the binary orbit where the ef- 
fective Roche lobes of the component stars are smallest 
(iLavton et al.lll998l:lSeDinskv et al.ll2007aHChurch et al.l 
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120091 : iRegos et al.ll2005| ) . Until recently, studies of mass 
transfer in eccentric binaries mainly focused on smoothed 
particle hydrodynamics calculations of the mass transfer 
stream and much less on the effects of the mass transfer 
on secular binary evolution. 

To the best of our knowledge. iSepinskv et al.l ()2007bL 
l2009f l were the flrst to study in detail the secular evo- 
lution of the orbital semi-major axis and eccentricity 
of eccentric mass-transferring binaries. The authors 
adopted a perturbati on method initially outlined by 
iHadjidemetrioul ()1969[ ) and found that, depending on the 
binary properties at the onset of mass transfer, the or- 
bital semi-major axis and eccentricity can either increase 
or decrease at a rate linearly proportional the mass trans- 
fer rate at periastron. Contrary to common assumptions, 
tides therefore do not always rapidly circularize binary 
o rbits as a star e v olves towards filling its Roche lobe. 

ISepinskv et al.l ()2007bL l2009f ) calculated the orbital 
evolution due to mass transfer in eccentric binaries as- 



suming mass transfer leads to the formation of an accre- 
tion disk. Tidal interactions between the disk and the 
donor star then allow at least some fraction of the an- 
gular momentum carried by the transferred matter to be 
returned to the orbit. However, if the size of the accretor 
is a considerable fraction of the size of the orbit, mat- 
ter lost by the donor star through the inner Lagrangian 
point may impact the accretor's surface directly rather 
than form an accretion disk. 
In this pap er, we complement the work of 

hereafter Paper I) and 

by investigating the orbital 

overflow in eccentric binaries 

the accretor. For this 
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evolution due to mass 
in which no disk forms around 
purpose, we calculate the trajectories of mass elements 
escaping through the inner Lagrangian point of the 
donor star in the ballistic limit, and distinguish between 
three possible cases: (1) the ejected matter falls back 
onto the donor star, (2) the ejected matter impacts the 
surface of the accretor, and (3) the ejected matter orbits 
the accretor for at least one orbital revolution of the 
binary. In the first two cases, the evolution of the orbital 
elements is obtained directly from the ballistic trajectory 
calculations. In the third case, the evolution of the 
orbital elements depends on the interactions between 
mass elements ejected during successive periastron 
passages of the donor star. Since such interactions are 



expected to violate our assumption of particles moving 
in the ballistic limit, we do not consider them here. 

The present paper is organized as follows. In §2, we 
outline the basic assumptions and methodology used to 
study ballistic particle trajectories and the effects of mass 
overflow on the orbital evolution of eccentric binaries. In 
§ 3, we examine the fate of mass lost by the donor star 
as a function of the initial binary properties and identify 
ranges of initial parameters leading to fallback, direct im- 
pact accretion, and possible disk formation. Evolution- 
ary sequences and timescales of orbital evolution due to 
mass overflow in eccentric binaries are presented in § 4. 
The final section, § 5, is devoted to concluding remarks. 

2. MASS OVERFLOW IN ECCENTRIC BINARIES 
2.1. Basic assumptions 

We consider a close binary system of stars with masses 
Ml and M2, and spherical radii TZi and TZ2- We de- 
fine the binary mass ratio as g = M1/M2, and assume 
the initial orbit to be Keplerian with semi-major axis a, 
eccentricity e, and orbital period Porb- The periastron 
distance of the two stars is then rp — a(l — e). 

The binary component stars are assumed to rotate 
uniformly with rotational angular velocities Qi and f22 
around axes perpendicular to the orbital plane and in the 
same sense as the orbital motion. For convenience, we 
express the rotational angular velocities in units of the 
orbital angular velocity Qp at periastron: 



ill — /i ^^p^ ^^2 — .72 ^^p 
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Hence, /i and /2 are a dimensionless parameterization of 
the rotation rates of stars 1 and 2, respectively, in terms 

oi\np\. 

For stars in a circular synchronous binary orbit, it is 
commonly assumed that the surface of the star will con- 
form to the cquipotcntial surface which encloses the same 
volume as the spherical radius of the star. The maxi- 
mum size of the star is then determined by the size of 
the Roche lobe, which is defined by the equipotential 
surface passing through the Li point. For eccentric bi- 
naries, we follow ISepinskv et al.l (2007a, b) and construct 
effective equipotential surfaces using a reference frame 
corotating with the star such tha t the envelope is in- 
stantaneously qu asi-static (see also lPlave3ll958t iLimbeii 
119631: IAvnil ll976D. For an analys i s of th e validity of this 
assumption, see ISepinskv et al.l (|2007a[ ). In this refer- 
ence frame we can locate the effective Li point at the 
minimum of the potential along the line connecting the 
mass centers of the two stars Q. We can then identify the 
effective Roche lobe as the equipotential surface passing 
through this effective Li point. 

At some point in time, one of the stars is assumed to 
fill its effective Roche lobe due to the expansion of the 

^ We note that viscosity in the stellar atmosphere will cause 
a tidal lag or lead in the position of Li with respect to line con- 
necting the centers of the two stars, to a degree dependent on the 
asynchronisni of the star. We intend to include tidal lags and leads 
in future investigations. 



star and/or the contraction of the orbit. For brevity, we 
refer to the effective Roche lobe-filling star as star 1 and 
to its companion as star 2. Since the stars are in an 
eccentric orbit, the size of star I's effective Roche lobe 
is a function o f time and is smallest when the stars are 
at periastron (jSepinskv et al.l [2007al hereafter SWK). 
Roche lobe overfiow is therefore expected to take place 
first at periastron and then re-occur at each subsequent 
periastron passage. 

To model the episodic mass loss at successive perias- 
tron passages, we assume a small amount of mass M3 
(M3 <C Ml, M2) to be ejected from star 1 through its in- 
ner Lagrangian point Li at each periastron passage. As 
in Paper I, we assume the mass overfiow rate Mi from 
the star to be a Dirac delta function of amplitude Mq < 
centered on the periastron of the binary orbit: 



Ml = Mo 5{v) 



(3) 



where v is the true anomaly with i' = at the perias- 
tron, and a normalization factor of 2 tt has been absorbed 
in the constant value for Mq. The mass of the ejected 
particle is then related to the mass overflow rate Mq by 



(1 - ef^ P, 
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M3= / {-Mi)dt = - 

-Po.b/2 {l + ey'^ 27r 



where we made use of the relation 



dt = 



(1 



e'f" Porb 



(1 + ecosi^) 



27r 



dv. 



Mo, (4) 
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In figure [ll we show the orbit for a system with ini- 
tial eccentricity e^ = 0.3, initial donor rotation rate 
fi^i ~ 0.8; and initial mass ratio g^ = 2.0 in a reference 
frame where star 2 is stationary. The ejected particle 
undergoes self-accretion, and its trajectory is shown by 
the dotted line. The orbit of the donor star is eccen- 
tric and shown by the short-dashed line. The position of 
star 2 is indicated by the point, and the dashed-dotted 
line shows its spherical radius. The donor star is shown 
at three positions (from right to left, respectively): peri- 
astron, particle impact, and apastron, and moves in the 
counter-clockwise direction. The spherical radius of the 
donor star, TZi, is shown by the long-dashed line and is a 
constant radius for each of the positions shown. The ef- 
fective Roche lobe of star 1 is also drawn at each position. 
It is clear that the effective Roche lobe radius of star 1 
grows relative to its initial size as it moves away from 
periastron. We also indicate the true anomaly, v, which 
measures the angle of the line connecting the mass cen- 
ters of the stars in a Keplerian orbit measured counter- 
clockwise from periastron. 

Furthermore, we consider the particle to impact the 
surface of a star when it falls within the spherical radius 
of that star. This surface does not lie at the Roche lobe, 
as can be seen in Figure [T] Due to the computational 
complexity and time- variable shape of the effective Roche 
lobe, calculations of this kind were not considered here. If 
this calculation were included, it is clear that the particle 
would strike the surface of the effective Roche lobe earlier 
in the orbit. In fact, for many of the cases studied, the 
ejected particle never moves further from the donor star 
than the position of the Li point at periastron. Thus, for 
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Fig. 1 . — A schematic representation of the orbit of the donor star (short-dashed Hne) in a reference frame fixed to star 2. In this diagram, 
we show the constant radius of star 2 (dashed-dottcd Unes), the constant radius of the donor (long-dashed hne), the time-dependent effective 
Roche lobe (solid line), and the path of the ejected particle (dotted line). The donor star is shown at three different positions in the orbit: 
periastron, particle impact, and apastron. The true anomaly, v, is the angle of the line connecting the mass centers of the stars measured 
with respect to periastron, and is indicated above. For this orbit (g^ = 2.0, fi^i = 0.8, e^ = 0.3), the ejected particle impacts the donor star 
at u ct 1.21 radians. 



these cases, if we assumed impact occurred at the surface 
of the effective Roche lobe, no mass loss would occur at 
the Li point, and thus we would not expect any orbital 
evolution due to this effect. As such, it may be possible 
for non-synchronous, eccentric systems to effectively halt 
mass loss through the Li point. We will investigate this 
effect in a future study. 

2.2. Ballistic particle trajectories 

Our aim in this paper is to investigate the evolution of 
the orbital elements due to mass overflow in eccentric bi- 
naries. For this purpose, we treat the binary components 
as rigid spheres of uniform density. The matter ejected 
by the donor star at each periastron passage is treated 
as a point mass. To describe the motion of the stars 
and the ejected mass elements, we introduce an inertial 
frame of reference OXYZ with origin O at the initial 
center of mass of the binary prior to any mass overflow, 
and with the XY-plane coinciding with the orbital plane. 
The orientation of the axes is chosen such that the posi- 
tive direction of the X-axis coincides with the line from 
the binary mass center to the position of the donor star 
at the onset of mass overflow, and the positive Z-axis is 
parallel to the direction of the orbital angular momentum 
vector at the onset of mass overflow. 

To determine whether a star undergoes mass over- 
flow, its volume-equivalent Roche lobe radius is calcu- 



lated accounting for the effects of eccentricity and non- 
synchronous rotation, as described in detail in SWK. 
Since we here do not consider the detailed structure of 
the donor star, we assume the matter ejected at peri- 
astron to have an initial velocity V3 equal to the vector 
sum of the star's rotational velocity at Li and its orbital 
velocity at periastron: 



Vs^^ni -xrr +V^i. 



(6) 



Here, r^^ is the vector from the center of mass of star 1 

to the Li point, and Vi the orbital velocity of star 1 at 
periastron right before the ejection of the matter. We 
do not consider any contribution to V3 from the thermal 
speed of the mass e lements in the donor's atmosphere 
(jLubow fc Shul II975D . A generalization of our calcula- 
tions accounting for the thermal speed is straightforward, 
but adds extra dimensions to the parameter space ex- 
plored in this paper. We therefore defer this extension 
to a future investigation. 

After the ejection of a mass element from the Li point, 
the system evolves as a three-body system according to 
the equations of motion 



df^ 



— E 



GM. 



rkj 



3 = 1 

J9^k 



r/cj-r Fkj] 



(fc = l,2,3), (7) 



where G is the Newtonian gravitational eonstant, i?i and 
i?2 are the position vectors of star 1 and star 2 with re- 
spect to the inertial frame of reference, i?3 is the position 
vector of the ejected mass element with respect to the in- 
ertial frame, and Vkj = Rk — Rj- 

For numerical convenience, we express equations ([7]) 
and the initial condition given by equation (jS]) in terms 

dimensionless quantities R*, M*, t* , V*, and (i*p defined 



R*=Rj/rp,i, j^ 1,2,3 
M* = M,/Mt, j = 1,2,3 

J = 1,2,3 




i Ip — i Ip 



(8) 
(9) 

(10) 
(11) 
(12) 



where rp- is the initial periastron distance between the 
binary components at the onset of mass overflow, Mt is 
the total system mass, and the subscripts j = 1, 2, 3 refer 
to star 1, star 2, and the ejected particle, respectively. In 
equation ^, we also express the position of the Li point 
in units of the initial periastron distance as 



(13) 



The quantity r£ depends only on the binary mass ratio, 
the orbital eccentricity, the true anomaly, and the rota- 
tion rate of star 1 in units of the orbital angular velocity 
at periastron (See Equation [A15] of Paper 1). 

In terms of these dimensionless quantities, the initial 
positions and velocities of the stars and ejected particle 
are 



o* — / V* v^* '7* \ / 
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{Xl,,Y*^„ Zl,) = (-T^, 0, 0), (15) 

(16) 
(17) 
(18) 
(19) 
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where the subscript "i" indicates that the quantity is to 
be considered at the time the particle is ejected, recalling 
that the particle is ejected at periastron. Furthermore, 
equation [7] becomes 
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(fc = 1,2,3). (20) 



In this form, both the initial conditions and the equations 
of motion are independent of the initial size of the orbit. 



We use an 8th order Runge Kutta ordinary differential 
equation solver (jGalassi et al.l [20061 ) to integrate equa- 
tions pO|) forward in time, starting from the first peri- 
astron passage at which matter is ejected by the donor. 
The three-body integration ends when (1) the ejected 
matter impacts star 2 ("direct impact accretion"); (2) the 
ejected matter falls back onto star 1 ("self-accretion"); or 
(3) the donor star has completed a full orbital revolution 
and a new matter element is ejected at the periastron of 
the binary orbit. After the accretion of the ejected mat- 
ter onto star 1 or 2 in cases (1) and (2), equations (O are 
integrated as a two-body problem by setting A/3 = 0. 

In binaries with eccentric orbits, fallback onto the 
donor star (case 2) is facilitated by the increasing size 
of the effective Roche Lobe as the distance between the 
stars increases following periastron passage (See SWK). 
The actual outcome of the mass ejection, though, also 
depends on the initial velocity of the ejected matter. In 
the case where no accretion takes place within one orbital 
revolution (case 3), the system evolves as a four-body 
system until one of the ejected matter elements is ac- 
creted by star 1 or star 2. or until another mass element 
is ejected at the next periastron passage. The successive 
ejection of mass elements during subsequent periastron 
passages, possibly causing the mass transfer stream to 
interact with itself, will likely lead to the formation of 
an accretion disk around star 2. A detailed study of ac- 
cretion disk formation will not be investigated here. We 
therefore restrict ourselves to identifying the regions of 
parameter space where this outcome is expected. 

Some example trajectories of particles ejected from the 
Li point of star 1 at the periastron of the initial binary 
orbit are shown in Figure [2l for binaries with different 
initial mass ratios, donor rotation rates, and orbital ec- 
centricities. The initial positions of the stars and the 
ejected matter are marked by asterisks, and the X* and 
Y* coordinates are expressed in units of the initial pe- 
riastron distance. The motion of the stars is counter- 
clockwise about the origin, and the numbered points on 
the orbits show the position of that object at different 
equally spaced values of the true anomaly. In all cases 
the radius of the donor star is set equal to the volume- 
equivalent radius of its effective Roche lobe at the perias- 
tron of the initial binary orbit (see SWK), and the mass 
of the ejected matter is determined by means of equa- 
tion (g]). We note that as long as M3 < Mi,M2, the 
ballistic trajectories are independent of the actual value 
of the mass of the ejected matter. For the presented 
trajectories, star 2 is assumed to be a point mass. 

We note that the shape of the particle trajectories, and 
hence the outcome of the mass overflow, is independent 
of the size of the orbit. This is because the equations 
of motion and the initial positions and velocities are all 
independent of the size of the orbit (see equations [l4] - 
ED]). 

In panels (a), (c), and (e), the ejected particles fall 
back onto the donor star in less than half an orbit, while 
in panels (b), (d), and (f) the particles end up orbiting 
the donor's companion. Since the companion is assumed 
to be a point mass with zero radius, none of the dis- 
played trajectories lead to direct impact accretion. In 
panel (d), direct impact would occur if the companion 
had a radius greater than 1.5% of the initial periastron 
distance. Comparison of panels (a) and (b) shows that 
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Fig. 2. — Example trajectories of particles ejected from the Li point in binaries with different initial mass ratios {qi — A/i^i/A/2,i), donor rotation 
rates (/i,i), and orbital eccentricities (ci). The X* and Y* coordinates are Cartesian coordinates in the orbital plane expressed in units of the 
initial periastron separation between the two stars and with origin at the center of mass of the binary system. The solid lines represent the path 
of the ejected particle, the dashed lines the path of the donor star (star 1), and the dotted lines the path of the companion star (star 2), which 
is assumed to be a point mass. The motion of the stars and particles is counter-clockwise about the origin, starting from the asterisks indicating 
the initial position of the component stars and the ejected matter at the onset of Roche lobe overflow. The numbered points show the position of 
the stars and the ejected particle at seven different equally spaced values of the true anomaly. The orbits of the stars are drawn for one full orbit, 
regardless of whether the particle has accreted or not. In panels (a), (c), and (e), the particle paths are shown up until the point where the particle 
is accreted back onto the donor star. In panels (b), (d), and (f), the particles undergo neither direct impact nor self-accretion. The paths of the 
particles are therefore shown for one complete binary orbit. 
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small changes in the orbital eccentricity can significantly 
affect the type and outcome of the particle trajectories. 
In this particular case, the slightly smaller initial veloc- 
ity of the ejected particle due to the lower eccentricity 
orbit in panel (b) is enough to cause the particle to fall 
towards star 2 rather than fall back onto star 1. Com- 
parison of panels (c) and (d) as well as panels (e) and (f ) 
shows that similar transitions in the behavior of the par- 
ticle trajectories can be obtained by changing the initial 
binary mass ratio or the initial donor rotation rate. 

2.3. Spin evolution 

Since the initial velocity of the matter ejected at pe- 
riastron depends on the rotation of the donor star, it is 
important to keep track of changes in the star's rota- 
tional angular velocity due to mass loss and accretion. 
To calculate the change in the rotational angular veloc- 
ity of the donor star due to the loss of an element of mass 
M3, we note that the total angular momentum Cj in the 
system right before the ejection of matter at periastron 
is given by 

+MijRi,i X Vij + AhjRij X V2,i, (21) 

where /i and I2 are the moments of inertia, i?i and R2 
the position vectors, and Vi and V2 the velocities of star 1 
and star 2, respectively, measured with respect to the in- 
ertial frame of reference. The subscript "I" denotes that 
the quantities are to be considered right before the ejec- 
tion of matter at periastron. Similarly, the total angular 
momentum Cp in the system right after the ejection of 
matter at periastron is given by 

Lp = ^l,F^l,_F + l2,F^2,F 

+ Mi,fRi,F X Vi,F + M2,fR2,F X V2,F 

+M3R3 X V3, (22) 

where R3 and V3 arc the position and velocity vectors of 
the ejected mass element at the time of ejection measured 
with respect to the inertial frame of reference. The sub- 
script "F" denotes that the quantities are to be consid- 
ered right after the ejection of matter at periastron. For 
the uniform density spheres considered here, the moment 
of inertia is / = OAMTZ^. For a more centrally condensed 
object the moment of inertia will decrease, reducing the 
contribution of the spin to the total angular momentum. 
Since the ejection of matter by star 1 does not instanta- 
neously affect the companion star, all quantities pertain- 
ing to star 2 are identical right before and right after the 
instantaneous ejection of a mass element. Conservation 
of total angular momentum therefore yields 

1 



(23) 



^i,F = -T- h.i^ij + MijRi.i X Vi.i 

J^I.F ^ 

-Mi^fRi.F X Vi^F - M3R3 X V3 



In this equation, all quantities on the right-hand side 
are readily obtained from the initial conditions and the 
adopted mass overflow rate, except for Ri^f and Vi^f- 
To determine Ri^f we assume M3 is ejected perfectly in- 



clastically, i.e., with no loss of linear momentund. Fur- 
thermore, since the ejection is instantaneous, the center 
of mass of star 1 before the ejection is identical to the 
center of mass of the system consisting of the particle 
and star 1 after the ejection. Thus, the final position of 
star 1 is given by 



^1, 



1 



Ml, 



MijRij - M3R. 



(24) 



The velocity Vi^f is determined from the conservation of 
linear momentum: 



ViJ 



1 



M- 



i,F 



Mi.iVij - M3V. 



(25) 



To calculate the change in the rotational angular ve- 
locity due to the fallback of the ejected matter onto the 
donor star, it suffices to replace M3 by — M3, and R3 and 
V3 by the position and velocity vector of the ejected mass 
clement at the time of accretion in equations (|23p - (P5|) . 

The change in the rotational angular velocity of star 2 
due to direct impact accretion is derived similarly: 



n 



2,F- 



1 
I2,F 



h.I^2J + M21R2.I X V2.1 



VM3R3 X Fg - A/2,Fi?2,F X V2A (26) 



with 



and 



R2.F = 



V^2,F = 



M2.F 



A/2,F 



M2JR2J + M3R3 



M2JV2J + M3V3 



(27) 



(28) 



Here, the subscripts "I" and "F" denote the quantities 
are to be considered right before and right after the ac- 
cretion of matter, respectively, and R3 and V3 are the 
position and velocity vector of the ejected mass element 
at the time of accretion. 

2.4. Calculation of the Orbital Elements 

Due to the gravitational interactions between the 
ejected matter and the component stars of the binary 
system, the orbital motion of the stars is no longer purely 
Keplerian. In particular, the semi-major axis and eccen- 
tricity are not constant in time. We calculate the instan- 
taneous value of the osculating orbital elements of the 
binary system at each timcstep using (ptcrnc 19601 ) 



= 1- 



G{Mi+M2)a 



1 



vr. 



(29) 



(30) 



a |f,ei| G(A/i + M2)' 

^ We note that, in a perfeetly inelastic collision energy is not 
conserved. We assume the existence of some mechanism, either 
internal or external, which allows the ejected particle to become 
physically separate from the donor star. Similarly for accretion, 
some process must exist to allow the particle and accretor to phys- 
ically merge. In either case, the change in the total energy of the 
system is likely to happen at the expense of the thermal energy of 
the component s tars. For further discussion of the change in the 
total energy, see i|4.2l 
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where fiei and v^ei are, respectively, the position and ve- 
locity vectors of star 2 with respect to star 1. The os- 
culating orbital elements represent the value of a and e 
that the orbit would have if the pcrturbative forces were 
removed at that instant. To keep track of changes in the 
position of the periastron due to mass overflow, we also 
calculate the instantaneous true anomaly of the binary 
system from 



1000 



1 

COS V ~ — 
e 



(1-^^) ; 



rrel 



(31) 



3. PARAMETER SPACE FOR DIRECT IMPACT AND 
SELF-ACCRETION 

In order to investigate the parameter space leading to 
different outcomes of mass overflow (direct impact ac- 
cretion, self-accretion, or possible disk formation) in bi- 
naries with eccentric orbits, we performed an extensive 
set of ballistic orbit calculations for binaries with differ- 
ent initial mass ratios g,, initial donor rotation rates /i_i, 
and initial orbital eccentricities e^ at the onset of mass 
overflow. As in ^2.2[ we set the donor's radius equal to 
the volume equivalent radius of its effective Roche lobe 
at the periastron of the initial binary orbit. We recall 
that as long as M3 ^ M\ , M2 the problem is fully de- 
termined by Qi, /i_i, and e^, and is independent of the 

actual values of the mass overflow rate Mq , donor mass 
Ml, donor radius, TZi, and initial orbital semi- major axis 
tti. Star 2 is again treated as a point mass. 

For each initial orbital eccentricity and donor rotation 
rate, we find a critical mass ratio, qait, such that for 
Qi > 9crit the system undergoes self-accretion, while for 
Qi < 9critj the system may undergo direct impact accre- 
tion or form an accretion disk, depending on the radius 
of star 2. The critical mass ratio is shown in Figure [3] as 
a function of the donor's initial rotation rate for initial 
orbital eccentricities ranging from 0.01 to 0.2. In this 
eccentricity range, gcrit decreases with increasing initial 
orbital eccentricity, increasing the parameter space lead- 
ing to sclf-accrction. 

For initial eccentricities e^ > 0.2, mass overflow al- 
ways leads to self-accretion if the initial binary mass ratio 
Qi > 0.006. In high-eccentricity systems direct impact ac- 
cretion therefore only occurs for very low mass ratios. For 
initial donor rotation rates in the range 1.0 < fi^i < 1.2, 
mass ejection at periastron results in self-accretion for all 
initial mass ratios lO^'^ < qi < 10^, independent of the 
initial orbital eccentricity. For the purpose of this paper, 
we did not thoroughly explore the parameter space for 
fi^i > 1.2, though calculation of some sample trajectories 
of mass ejected from highly super-synchronously rotat- 
ing donor stars indicates the matter may start orbiting 
the donor star before leading to self-accretion. 

When Qi < (7crit, depending on the accretor's radius 
72-2 the ejected matter can directly impact the compan- 
ion star or go into orbit around it, possibly forming an 
accretion disk. While studying disk formation is beyond 
the capabilities of our ballistic trajectory code, we can 
investigate when direct impact rather than disk forma- 
tion is expected to occur based on the ballistic trajecto- 
ries of the ejected mass elements. For this purpose, we 
treat star 2 as a point mass and determine the minimum 
distance rmin between the ejected matter and star 2 dur- 
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Fig. 3. — The critical mass ratio 5crit above which systems un- 
dergo self-accretion as a function of the donor's initial rotation rate 
/l i, for initial orbital eccentricities ranging from 0.01 to 0.2. For 
Qi > Qcrit) systems undergo self-accretion, while for qi < g^rit sys- 
tems undergo direct impact accretion or form an accretion disk, 
depending on the radius of the accretor. We note that for low ini- 
tial eccentricities, e^ < 0.05, the delta-function approximation of 
the mass overflow rate given by equation Jsjl likely breaks down, 
so that at low eccentricities disk formation or direct impact may 
occur for a larger region of parameter space than presented here. 

ing a single orbit for ballistic trajectories not leading to 
sclf-accrction. Hence, if 7^2 > ''min; direct impact oc- 
curs within one binary orbital period. If 7^2 < ''min, the 
ejected matter may still impact star 2 after one orbital 
period, but it may also orbit star 2 until it interacts with 
other ejected mass elements and forms an accretion disk. 

The minimum distance r^in/fp,i is shown in Figure |4] 
as a function of the logarithm of the initial binary mass 
ratio Qi. For each initial mass ratio, the grey dots rep- 
resent the minimum distance associated with different 
initial donor rotation rates (0.01 < fi^i < 1.0) and ini- 
tial orbital eccentricities (0.01 < e^ < 0.2). The main 
trend of rmin/fp.i as a function of qi is to decrease with 
increasing values of qi. This is partly because the Roche 
lobe of star 2 is smaller for large mass ratios than for 
small mass ratios, or, equivalently, the initial position of 
the ejected matter, the Li point of star 1. is relatively 
closer to star 2 for large mass ratios than for small mass 
ratios (see SWK). 

The lower and upper edges of the grey dotted region 
in Figure HI can be approximated by 



^min, lower 



rp. 



exp [-0.16(logio q^ + 4.8)2 - 0.13] , (32) 



and 



' min, upper 



rp,i 



= exp[-0.1(logiog. + 4.8)2] 
- 0.01 logiog, + 0.05. 



(33) 

Hence, if 7?.2 > rmin, upper, the system always undergoes 
direct impact accretion within one orbital period. If 
7?-2 < J'min, lower, thc cjccted matter always avoids im- 
pacting the donor's companion star within one orbital 
period, so that an accretion disk may be formed. If 



.lower 



< 112 < rn 
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Fig. 4. — The minimum distance rmin/rp i from star 2 (assumed 
to be a point mass) reached by matter ejected from the Li point 
of the donor star at the periastron of the binary orbit as a function 
of the logarithm of the initial binary mass ratio qi. The distance 
is measured in the orbital plane and is normalized to the initial 
periastron separation rpi between the two stars. For a given mass 
ratio, the grey dots represent minimum distances for different ini- 
tial donor rotation rates /i ^ ranging from 0.01 to 1.0, and initial 
orbital eccentricities e; ranging from 0.01 to 0.2. The solid lines 
represent the approximate upper and lower limi ts o n the minimum 
distance as a function of qi g iven by equations II32II and l|33| l . The 
black circles are the data of ILubow fc Shul II1975I ) for the closest 
approach of a ballistic particle in a circular, sync hronous binary 
system , and the dashed line is the fit t o the data of lLubow &: Shul 
II1975I ) given bv lNelemans et al.l 1 120011 ). Both the fit and the data 
shown are adjusted to account for the different definitions of the 
mass ratio given in those papers. 

overflow depends on the initial donor rotation rate and 
orbital eccentricity. 

For comparison, the black circles of Figure |4] show 
the distance of closest approach for a particle in a 
circular, synchronou s binary system as calculated by 
ILubow fc Shul (|1975[ ) (their Wmin), where we have taken 
into account the difference between their and our def- 
inition of the mass ratio. We note that this data is 
consistent with our calculation of the distance of clos- 
est approach, and that the inclusion of non-synchronous 
rotation and/or eccentricity can either increase or de- 
crease the distance of closest approach. The dashed 
line through t h ese d ata points is the fit to the data o f 
ILubow fc Sliul p975l ) given by iNelemans et~an (J200lb . 
The fit is very accurate for gi < 2, but should not be 
extrapo lated to large q j . Th is does not affect the re- 
sults of INelemans et aLl (|2001l ). as they are interested in 
double white dwarf systems where the mass ratio of the 
system is generally less than unity. 

4. ORBITAL EVOLUTION DUE TO MASS OVERFLOW 

Using the framework described in ^ we can now inves- 
tigate the secular evolution of the orbital elements due to 
the effects of mass overflow in eccentric binary systems. 
Combining equations (f23| -([28 |) with equations (|29|) -([3T |) . 
it can be seen that the changes in the orbital elements, 
and thus the time scales of orbital evolution, depend on 
the binary component masses, the donor radius and ro- 
tation rate (through the initial particle velocity) , and the 
orbital elements (also through the initial particle veloc- 



ity). The outcome of the mass overflow (e.g. whether or 
not direct impact accretion occurs) also depends on the 
radius of the accretor. 

Once mass overflow starts, the radius of the donor star 
evolves due to the combined effects of mass loss and stel- 
lar evolution. In addition, the radius of the donor's ef- 
fective Roche lobe changes due the evolution of its ro- 
tational angular velocity, the orbital elements, and the 
binary mass ratio (SWK). In order for mass overflow 
in eccentric binaries to be stable and long-lived, the ra- 
dial evolution of the donor star must follow the evolu- 
tion of its effective Roche lobe at the periastron of the 
binary orbit. In this section, we explore the potential 
long-term evolution of effective Roche lobe overflowing 
binaries with eccentric orbits, assuming that after the 
onset of mass overflow the donor star's radius remains 
equal to the radius of its effective Roche lobe at each 
subsequent periastron passage. Whether the response of 
the star to mass loss renders this possible depends on its 
internal structure and on the binary properties. Since 
we here treat the stars as rigid spheres, we defer ques- 
tions on the longevity and the stability of mass overflow 
to future investigations. 

To illustrate the effects of mass overflow in eccentric bi- 
naries, we consider a binary star system consisting of two 
rigid spheres with radii representative of main-sequence 
stars. The first rigid sphere ("star 1") is assumed to 
have a radius representative of a star near the end of 
core-hydrogen burning when the main-sequence radius 
is maximal, while the second rigid sphere ("star 2") is 
assumed to be a zero-age main-sequence star. Such a bi- 
nary consisting of two non-coeval component stars can be 
formed through dynamical interactions in globular clus- 
ters. The adopted configuration allows us to illustrate 
the orbital evolution time scales of binaries with initial 
mass ratios qi < 1 as well as binaries with initial mass 
ratios qi > 1. In addition, a zero-age main-sequence 
accretor constitutes a big enough "target" to provide a 
non-negligible parameter space leading to direct impact 
accretion. A study of direct impact accretion in compact 
object binaries, double white dwarfs in particular, will 
be presented in a forthcoming paper. 

4.1. Orbital Evolution Timescales 

In Figures [5]-[51 we show the variations of the instanta- 
neous characteristic timescales r^ = \a/a\ for the evolu- 
tion of the orbital semi-major axis a and t^ = |e/e| for the 
evolution of the orbital eccentricity e, as contour plots in 
the {qi, ei)-plane at the start of mass overfiow. The con- 
sidered range of initial mass ratios is 0.01 < qi < 100. 
Each panel in the figures corresponds to different values 
of the initial donor mass Mi^i and the initial rotation rate 
/i.j. The mass overflow rate at periastron is taken to be 
A/q = — 10^^ Af© yr~^. Since the resulting timescales 
are inversely proportional to Mq, the results can easily 
be rescaled to different mass overflow rates. The radii of 
the two stars are calculated us ing the ana lytic formulae 
for stellar evolution derived bv lHurlev et al. (200 0i). 

To calculate the instantaneous rates of secular evolu- 
tion of the orbital elements, we follow the evolution of 
the orbital semi-major axis and eccentricity of eccentric, 
effective Roche lobe overflowing binaries for 100 consec- 
utive orbits. During this evolution, transitions may take 
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Fig. 5. — Timescales for the evolution of the semi-major axis a as a function of the initial mass ratio qi and initial orbital eccentricity e^ 
in systems undergoing direct impact accretion. The donor star is assumed to be near the end of the main sequence and to have an initial 
mass Mi^i = 1,3, or 10 Mq, and an initial rotation rate fi^i = 0.5,0.75, or 0.9. The accretor is assumed to be a zero-age main-sequence 

star with a mass determined by the initial mass ratio qi. The mass overflow rate Mo at periastron is — 10~^ Mq yr~^. The different shades 
of gray represent regions in the {qi,ei) parameter space with different characteristic timescales of orbital evolution. From the darkest to 
the lightest shade of gray, the timescales are: OGyr < \Ta\ < 1 Gyr, 1 Gyr < \Ta\ < 5 Gyr, 5 Gyr < \Ta\ < 10 Gyr, lOGyr < \Ta\ < 15Gyr, 
15Gyr < |ra|. The dotted black line near g^ = 1 divides the regions of the {qi, e^) parameter space where the semi-major axis increases as 
a function of time (to the left of the line) from the regions of the (qi, ei) parameter space where the semi- major axis decreases as a function 
of time (to the right of the line). The black shaded area corresponds to regions of the parameter space where the mode of mass overflow 
switches from direct impact to self-accretion or possible disk formation within the first 100 orbits after the onset of mass overflow. This 
area occupies only a small space between the the direct impact and self-accretion regimes. Unshaded regions correspond to regions of the 
parameter space where mass overflow either does not lead to direct impact accretion or the zero-age main-sequence accretor is so large that 
a contact binary rather than a semi-detached binary is formed. 



place from direct impact accretion to self-accretion or 
vice versa. Systems for which this happens within 100 
orbits from the start of mass overflow arc confined to a 
fairly restricted region of the parameter space and are 
shown in black in Figures [SHSl In addition, for small ini- 
tial mass ratios qi < 0.01 — 0.02. the radius of star 2 is 
so large relative to the size of the orbit that the binary 
becomes a contact binary rather than a semi-detached 
binary. These binaries are also left out of consideration 
in Figures [SHHl For binaries that do not transition be- 
tween direct impact and self-accretion and that do not 



form contact binaries, timescales of orbital evolution are 
calculated by averaging the changes in the orbital ele- 
ments over one orbital period. 

Figures [5] and [S] show the orbital evolution timescales 
for the semi-major axis and eccentricity for the case of 
direct impact accretion. The dotted line in the figures 
separates regions of the parameter space where a and 
e increase (to the left of the dotted line) from regions 
of the parameter space where a and e decrease (to the 
right of the dotted line). For small orbital eccentricities, 
the transition from increasing to decreasing orbital semi- 
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Fig. 6. — Same as Figure [5] but for the evolution of the orbital eccentricity e. 



major axis occurs near qi = 1, consistent with what is 
expected for conservative mass transfer in binaries with 
circular orbits (see, e.g., Paper I). As the eccentricity 
increases, the dividing hne between increasing and de- 
creasing semi- major axis moves toward smaller qi. For a 
given choice of binary parameters, the orbital eccentricity 
typically evolves on shorter timescales than the orbital 
semi-major axis. The eccentricity evolution furthermore 
tends to circularize the orbit on fairly short timescales 
for most of the explored parameter space. Eccentricity 
pumping is only found for low initial binary mass ratios 
(qi < 0.1), low initial orbital eccentricities (e^ < 0.05), 
and substantially subsynchronous rotation (fi^i ^ 0.75). 
Figures [7] and [51 show the orbital evolution timescales 
for the orbital semi-major axis and eccentricity in the 
case of self-accretion. Unlike the direct impact case, self- 
accretion always acts to decrease both the semi-major 
axis a and the orbital eccentricity e. The evolution of 
both a and e is furthermore fastest for larger binary mass 
ratios and lower orbital eccentricities. Comparison of 



Figures [7] and [8] with Figures [5] and |6] furthermore shows 
that the contours for self-accretion systems complement 
those for direct impact systems. 

4.2. Example Orbital Evolution Sequences 

In Figures [9] and [10] we show the long-term evolu- 
tion of orbital and stellar properties for a representa- 
tive sample of systems undergoing direct impact and 
self-accretion, respectively, for a mass overflow rate of 
Mq = — 10~^ Mq yr"-'^. Following the procedure outlined 
in [21 a particle is ejected from the donor star at each 
periastron passage of the binary, and its ballistic orbit 
is followed until direct impact or self-accretion occurs. 
Impact of the particle on either star is determined un- 
der the assumption that the star is spherically symmetric 
and uniform in density (/ = OAMTZ^). Hence, we do not 
account for deviations from spherical symmetry due to 
rotation and tidal interactions. 

Following the accretion of the particle, the system 
evolves as a two body system until it reaches the pe- 
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Fig. 7. — Timescales for the evolution of the semi-major axis a as a function of the initial mass ratio qi and initial orbital eccentricity 
a in systems undergoing self-accretion. The donor star is assumed to be near the end of the main sequence and to have an initial mass 
A^l.i = Ij 3, or 10 Mq, and an initial rotation rate /i_i = 0.5, 0.75, or 0.9. The binary companion is assumed to be a zero-age main-sequence 



star with a mass determined by the mass ratio qi. The mass overflow rate at periastron is Mq 



-lO-^Mgyr- 



The different shades 



of gray represent regions in the {qi,ei) parameter space with different characteristic timescales of orbital evolution. From the darkest to 
the lightest shade of gray, the timescales are: OGyr < |t(,| < 1 Gyr, 1 Gyr < \Ta\ < 5 Gyr, 5 Gyr < \Ta\ < lOGyr, lOGyr < |t(,| < 15 Gyr, 
15Gyr < |ra|. The timescales always represent a decrease of the semi-major axis. The black points at low eccentricity correspond to 
regions of the parameter space where the mode of mass overflow switches from self-accretion to direct impact accretion or possible disk 
formation within the first 100 orbits after the onset of mass overflow. Those points occupy only a small space between the the direct impact 
and self-accretion regimes. Unshaded regions correspond to regions of the parameter space where either mass overflow does not lead to 
sclf-accrction or the zero-age main-sequence accrctor is so large that a contact binary rather than a semi-detached binary is formed. 



riastron of the new orbit. At this point, the effective 
Roche lobe radius will also have changed compared to 
its value at the time of the previous mass ejection due 
to changes in the binary mass ratio, semi-major axis, ec- 
centricity, and donor rotation rate. Since we assume the 
radius of the donor star to remain equal to the radius 
of its effective Roche lobe at periastron during the life- 
time of the mass overflow phase, the donor rotation rate 
must be updated in accordance with conservation of spin 
angular momentum before the next mass ejection at pe- 
riastron. Because the donor rotation rate in turn affects 
the radius of the effective Roche lobe, a root-finding al- 



gorithm is used to calculate a new self-consistent radius 
and rotation rate for the donor star immediately before 
the ejection of matter at each consecutive periastron pas- 
sage. Throughout our numerical integrations, the total 
angular momentum of the system (defined as the spin 
angular momenta of stars 1 and 2 plus the orbital angu- 
lar momentum of the system), is conserved to within a 
factor of 10~^. 

In Figures [91 andfTUl the fractional changes of the donor 
radius TZi, total energy, i^tot, eccentricity, e, semi-major 
axis, a, spin orbital angular momentum of stars 1 and 
2, Jsi and J52, and the rotation rates of stars 1 and 2, 
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Fig. 8. — Same as Figure [7] but for the evolution of the orbital eccentricity e. 



TABLE 1 
Initial and final binary parameters for the evolutionary sequences shown in Figures [g] and [TU] 





Ml 


(A%) 


M2 (Mq) 


9 




a(RQ) 


e 




/i 




/2 




tflnal 


System 


Initial 


Final 


Initial 


Final 


Initial 


Final 


Initial Final 


Initial 


Final 


Initial 


Final 


Initial 


Final 


(Gyr) 


A 


1.00 


0.78 


20.00 


20.02 


0.05 


0.04 


9.31 11.65 


0.025 


0.000 


0.50 


0.20 


1.00 


1.50 


0.14 


B 


1.00 


0.81 


10.00 


10.19 


0.10 


0.08 


7.85 9.63 


0.050 


0.000 


0.50 


0.23 


1.00 


1.58 


0.13 


C 


1.00 


0.99 


0.25 


0.26 


4.00 


3.81 


3.24 3.12 


0.050 


0.000 


0.75 


0.75 


1.00 


2.00 


0.01 


D 


3.00 


2.38 


15.00 


15.62 


0.20 


0.15 


2.26 28.56 


0.050 


0.000 


0.75 


0.31 


1.00 


2.49 


0.41 


E 


10.00 


9.97 


1.00 


1.03 


10.00 


9.68 


19.78 19.32 


0.025 


0.000 


0.90 


0.91 


1.00 


2.57 


0.02 


F 


1.00 


1.00 


0.25 


0.25 


4.00 


4.00 


3.16 3.07 


0.050 


0.000 


0.50 


0.53 


1.00 


1.06 


0.01 


G 


3.00 


3.00 


5.00 


5.00 


0.60 


0.60 


2.71 2.41 


0.400 


0.234 


0.75 


0.78 


1.00 


1.39 


1.42 


H 


3.00 


3.00 


0.50 


0.50 


6.00 


6.00 


17.41 14.23 


0.400 


0.063 


0.90 


0.88 


1.00 


1.66 


0.32 


I 


10.00 


10.00 


100.0 


100.0 


0.10 


0.10 


266.55 253.65 


0.800 


0.789 


0.50 


0.58 


1.00 


1.02 


11.9 



Note. — Mi is the mass of the donor star, M2 the mass of the companion, q = M1/M2 the mass ratio, a the semi-major axis, e 
the orbital eccentricity, /i the rotational angular velocity of the donor star in units of the orbital angular velocity at periastron, /2 the 
rotational angular velocity of the donor star in units of the orbital angular velocity at periastron, and tfi„al the time at which the calculation 
was terminated. 
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Fig. 9. — Evolution of orbital parameters as a function of time for systems undergoing direct impact accretion. The y-axes show the 
fractional change in the parameters relative to their initial value. The initial and final values for some of the system parameters are shown in 
table [T] In (a), we show, from top to bottom, the radius of star 1, the total energy of the system (orbital plus the spin of each component), 
the eccentricity, and the semi-major axis. In (b), we show, from top to bottom, the spin angular momentum of star 2, the spin angular 
momentum of star 1, and the rotation rates of stars 2 and 1. 
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/i and /2, arc shown as a function of time for systems 
undergoing direct impact accretion and self-accretion, re- 
spectively. The evolutionary sequences were calculated 
for up to four mass transfer timescales Mi/\Mq\ (sys- 
tems G and I) or until the binary circularized (systems A 
through F), whichever occurred first. For system H, the 
calculation was terminated at time t — 0.32 Gyr when it 
ceased to be a self-accretion system. The initial and final 
parameters of the nine evolutionary sequences shown are 
summarized in Table [T] 

The direct impact systems shown in Figure IHl all circu- 
larize within a few Gyr or less after the start of mass over- 
flow. In the case of system A, the eccentricity initially 
grows, reaching a peak emax = 0.0268 at i = 0.034 Gyr, 
before circularizing at an increasingly rapid rate. The 
semi-major axes generally increase, except for systems C 
and E, as expected for binaries with mass ratios q > 1 
and almost circular orbits. Mass loss from the donor 
star and the growth of the donor's radius, which is re- 
quired to track the growing effective Roche lobe, further- 
more decrease the donor's rotation rate to significantly 
sub-synchronous values at the end of the evolutionary 
sequence for systems A, B, and D. For systems C and 
E, mass transfer to the less massive star 2 causes the 
semi-major axis to decrease, increasing the orbital an- 
gular velocity at periastron, and thereby effectively in- 
creasing the donor rotation rate. This in turn causes the 
effective Roche lobe radius of the donor to shrink. We 
note that in all five cases shown, the end of the evolution- 
ary sequence corresponds to the time at which the binary 
becomes circular, not to the cessation of mass transfer in 
the binary system. 

The self- accretion systems shown in Figure [TOl decrease 
their eccentricity as a result of mass overflow, though sys- 
tem I has an evolutionary timescale for the orbital eccen- 
tricity much longer than a Hubble time (see Figure |8]). 
The orbital semi-major axes always decrease, while the 
donor rotation rates can increase as well as decrease dur- 
ing the course of mass overflow. There is no change in the 
spin angular momentum of star 2 since there is no mass 
transferred. In this case, the change in /2 is due entirely 
to the change in the orbital angular velocity at perias- 
tron. In the case of system G, the donor initially spins 
down by a few per cent relative to the rotational angular 
velocity at periastron, but reverses its spin evolution af- 
ter about 100 Myr. At the end of the evolutionary phase 
shown (corresponding to four mass transfer time scales) , 
the donor is spinning super-synchronously at the peri- 
astron of the binary orbit. The donor rotation rate for 
System H slowly decreases until just past t « 0.3 Gyr, at 
which time it rapidly plunges to very low rotation rates. 
We note that the physical rotation rate of this donor star 
never actually decreases, and that the decrease in /i is 
caused by the orbital angular velocity at periastron in- 
creasing at a rate faster than the actual rotation rate 
of the donor star. The rapid decrease starts a few or- 
bits before the system transitions from a sclf-accrction 
system to a system in which the ejected matter impacts 
neither star 1 nor star 2 during the course of one binary 
orbit. Since our code is currently not equipped to deal 
with the evolution of systems like this, we terminated the 
calculation at this point. 

We note that at every step of the orbital evolution 



calculations, both the total linear and angular momenta 
of the system remain constant. However, mass overflow 
can exchange momentum between the spins of the bi- 
nary components and the orbit. We particularly find 
that, contrary to pas t assmirpti ons (for exam ple, see 
IVerbunt fc RappaportI [1988 : Marsh et all |2004[ ). direct 
impact accretion does not necessarily provide a sink o f 
orbital angular momentum (see also iMotl et ahl |2007() . 
Instead, the transferred matter contains both spin and 
orbital angular momentum from the donor star, part of 
which can be returned to the orbit upon accretion. This 
makes it possible for the orbital angular momentum to 
increase at the expense of the spin angular momentum 
of the donor. The impact of this result on the stability 
of mass overflow in close binaries such as double white 
dwarfs will be the subject of a forthcoming investigation. 
Lastly, we can sec from Figures [9] and [10] that the to- 
tal energy of the system can change significantly over 
the course of its lifetime. Due to our assumption that 
the particle is ejected and accreted perfectly inelastically 
(i.e., with no loss of momentum), it is impossible to con- 
serve the total system energy as well. Even in the case 
of self-accretion where the mass of star 1 is unchanged 
the ejected particle self-accretes with a different position 
and velocity than that with which is was ejected. For 
the systems tested here, the total energy can change by 
as much as 10% or more, and can either increase or de- 
crease depending on the system parameters. The energy 
added or subtracted from the system might be accompa- 
nied by a commensurate change in the thermal energy of 
the stars' envelopes. 

5. DISCUSSION AND CONCLUSIONS 

In this paper, we extended our previous work on mass 
transfer in eccentric binaries by using ballistic trajectory 
calculations to determine the parameter space for differ- 
ent outcomes of mass overffow. Assuming mass overflow 
takes place through Roche lobe overflow at the periastron 
of the binary orbit, we explored a broad parameter space 
to determine the conditions under which mass overflow 
leads to direct impact accretion onto the donor's compan- 
ion or to fallback onto the donor star ("self- accretion"). 
The resul ts presented in this pap er along with those pre- 
sented in lSepinskv et all (|2007allbl . [2009I) . provide a self- 
consistent picture of the orbital evolution of eccentric 
effective Roche lobe overflowing binary star systems as 
a function of their initial orbital parameters. Until now, 
such a picture has been lacking both in binary stellar 
evolution and binary population synthesis codes. 

We find that systems with a large initial eccentricity 
(ci > 0.2) or a donor initially rotating near synchronic- 
ity [fi.i ^ 0.8) are always expected to undergo self- 
accretion. Direct impact accretion or disk formation is 
expected to occur mainly for systems with low eccentrici- 
ties, low mass ratios, and substantially subsynchronously 
rotating donor stars. Furthermore, self-accretion is found 
to always decrease the orbital eccentricity, and can do so 
over timescales ranging from less than a Myr to more 
than a Gyr depending on the initial binary parameters 
and the chosen mass overflow rate. 

To illustrate the applicability of the presented formal- 
ism, we calculated the orbital evolution due to mass over- 
flow for eccentric binaries consisting of an evolved main 
sequence donor and a zero-age main sequence accretor as 
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a function of the initial binary parameters. For binaries 
undergoing direct impact accretion, mass transfer can 
increase as well decrease the initial orbital semi-major 
axis and eccentricity, while for binaries undergoing self- 
accretion, mass overflow always decreases both the ini- 
tial orbital semi-major axis and eccentricity. For a mass 
overflow rate of Mq ~ — 10~^ Mq yr"-'^, the time scales of 
orbital evolution can range from less than 1 Gyr to more 
than a Hubble time. 

In this exploratory study, both the donor star and the 
accretor were treated as rigid spheres. The results pre- 
sented therefore do not account for orbital evolution due 
to tidal torques or for the radial response of the donor 
star to mass loss. Our calculations nevertheless clearly 
illustrate that mass overflow leading to direct impact ac- 
cretion or self-accretion can both enhance or counteract 
the effects of tides in eccentric effective Roche lobe over- 
flowing binarie s . A similar conclusion was reached by 
iSepinskv et al.l ()2007bl I2009D for binaries in which mass 
transfer leads to the formation of an accretion disk. As 
one of the next steps in our investigation, we will imple- 
ment more realistic stellar models in our ballistic particle 
trajectory code and examine the effects of tides and stel- 
lar evolution on mass transfer in eccentric binaries. 

While not discussed in detail here, we also find that di- 



rect impact accretion does not necessarily deposit all of 
the accreted matter's angular momentum into the spin 
of the accretor. Instead, some of this angular momen- 
tum, which originated from both spin and orbital an- 
gular momentum of the donor star, is deposited into 
the accretor's spin while the remaining part is returned 
to the orbit. This is in stark contrast to commonly 
adopted assumptions in stability studies of mass trans- 
fer iri ultra-compact binar ies fe.g.. lVerbunt fc Rappaporq 
119881 iMarshet all 120041 ). The possibility of avoiding 
strong orbital angular momentum loss during direct im- 
pact accretion is particularly relevant in predicting ex- 
pected gravitational wave detection rates from mass- 
transferring double white dwarfs (AM CVn stars) by the 
Laser Interferometer Space Antenna (LISA). 
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